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We present a study of transport in multiple-band non-interacting Fermi metallic systems based on 
the Keldysh formalism, taking into account the effects of Berry curvature due to spin-orbit coupling. 
We apply this formalism to a Rashba 2DEG ferromagnet and calculate the anomalous Hall effect 
(AHE) and anisotropic magnetoresistance (AMR). The numerical calculations reproduce analytical 
results in the metallic regime revealing the crossover between the skew scattering mechanism dom- 
inating in the clean systems and intrinsic mechanism dominating in the moderately dirty systems. 
As we increase the disorder further, the AHE starts to diminish due to the spectral broadening of 
the quasiparticles. Although for certain parameters this reduction of the AHE can be approximated 
as Oxy ~ o-^x with (p varying around 1.6, this is found not to be true in general as axy can go 
through a change in sign as a function of disorder strength in some cases. The reduction region in 
which the quasiparticle approximation is meaningful is relatively narrow; therefore, a theory with a 
wider range of applicability is called for. By considering the higher order skew scattering processes, 
we resolve some discrepancies between the AHE results obtained by using the Keldysh, Kubo and 
Boltzmann approaches. We also show that similar higher order processes are important for the 
AMR when the nonvertex and vertex parts cancel each other. We calculate the AMR in anisotropic 
systems properly taking into account the anisotropy of the non-equilibrium distribution function. 
These calculations confirm recent findings on the unreliability of common approximations to the 
Boltzmann equation. 

PACS numbers: 72.15.Eb, 72. 20. Dp, 72. 20. My, 72.25.-b 



I. INTRODUCTION 

Recently, the interest in transport calculations in 
multiple-band systemaii^ has been rekindled in part due 
to the realization of diluted magnetic semiconductors 
(DMS) that have strong spin-orbit interactions, vari- 
able carrier densities, and ferromagnetic ordering. These 
properties imply the existence of the anomalous Hall 
effect (AHE)^ and the anisotropic magnetoresistance 
(AMR)i^ Even though the mechanisms of the AHE and 
the AMR are different, they both have a similar descrip- 
tion based on the multiple-band transport theory. In this 
paper, we formulate a relatively simple framework for do- 
ing such transport calculations. 

The AHE is usually described in terms of the anoma- 
lous Hall resistivity pxy that measures the transverse 
voltage with respect to the transport direction and de- 
pends on the spontaneous magnetization M along the z 
direction. Theoretical studies of the AHE have a long his- 
tory beginning with the work of Karplus and Luttinger.— 
A number of papers on the AHE also appeared not so 
long ago^'^'^ d°i^^i^^ after the interpretation of the AHE 
based on the Berry phasei^ was proposed. Neverthe- 
less, theoretical description of the AHE is far from be- 
ing complete and it often involves cumbersome calcu- 
lations without transparent interpretationsj^ The diffi- 
culties appear due to the necessity to consider the off- 
diagonal elements in Bloch band indices (the interband 
coherences induced by charge currents). There is a gen- 



eral trend to focus on particular simple models in order to 
overcome the common mistakes that are made in treat- 
ing the AHE. A number of recent publications concen- 
trate on the simpler but non-trivial Rashba 2D electron 
system^ i-^^d^d'^i-^^d^i^°i^^i^^i^^ yet arriving at contradic- 
tory predictions. Most of the disagreements have been 
finally resolved22i>^>2i with some being addressed in this 
paper. 

In calculating the AHE for a given material, the usual 
approximations performed to leading order in h/rsF can 
fail, where r is the scattering time and is the Fermi 
energy. The semiclassical description of the Hall conduc- 
tivity within the usual Boltzmann equation leads to an 
AHE contribution due to the scattering asymmetry in the 
collision term usually labeled as skew scatteringi^ Other 
terms, arising from subtle issues dealing with interband 
coherence during the coUision and acceleration by the 
electric field between collisions, are usually introduced 
by hand through the so called anomalous velocity2£ and 
side-jumpi^ This approach however, is non-systematic 
and prone to errors from missing terms and wrong in- 
terpretations, e.g. such as giving physical meaning to 
gauge dependent quantities. A more systematic way to 
derive the correct semiclassical equations is through the 
Keldysh formalism in which these interband coherences 
effects are taken into account automatically-^i'^i 

The system under consideration also allows us to study 
the diagonal resistance as a function of the direction of 
the magnetization. The change in the resistance as a 
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function of the magnetization direction relative to the 
current or crystallographic direction is called the AMR 
effect. The microscopic origin of the AMR in transition 
metal ferromagnets is still elusive22ii2£i2£i2i and detailed 
calculations require consideration of complicated band 
structures . ^^1^^ A relatively simple host band structure in 
the DMS ferromagnets provides a possibility for perform- 
ing detailed microscopic calculations based on simple 
physical modelsi^ However, the relaxation time approx- 
imation used in such calculations is not always reliable 
since it does not fully take into account the anisotropics 
of the system!^ The Kubo formula approach has been 
applied to the AMR calculations in Rashba systems and 
it has revealed the cancellation of the nonvertex and ver- 
tex parts,— similar to the spin Hall effect (SHE) and the 
AHE. 

In this paper, we apply the Keldysh formalism for 
transport calculations in multiple-band non-interacting 
Fermi systems. This treatment simultaneously takes into 
account the Berry curvature effects (interband coher- 
ences) and scattering, allowing us to immediately account 
for such physical effects as side-jump scattering and skew 
scattering within the same footing. We calculate the 
AHE analytically and numerically for the Rashba model 
and find in agreement with Onoda et ali^^i^ three distinct 
regimes: the skew scattering regime, the disorder inde- 
pendent regime, and the dirty regime in which, although 
the basis of theory is not as well established, a distinct 
rapid reduction of the AHE is observed as the conductiv- 
ity o-xx diminishes. Even though almost all ferromagnetic 
systems are three dimensional, the findings of this simple 
2D model has been linked to higher dimensional systems 
arguing that most likely the major contributions to the 
AHE come from the band anti-crossing regions^ similar 
to one observed in the Rashba model. 

We further analyze the scaling found in the dirty 
regime^i^ in which the AHE seems to diminish in a man- 
ner that can be approximated as: axy ~ cr'^x with ip being 
close to 1.6. Some experimental results claim to confirm 
such scaling j^LSSii^iiflidi however, treatment of some of 
these experimental results has to be done with extra care 
as the region of interest is often restricted to less than a 
single decade, the materials have strong mangetoresis- 
tances and in-plane anisotropics associated with them, 
and most of the data associated with the zero field cal- 
culation is in fact at very high magnetic fields. Although 
our numerical results confirm this scaling, it is found to 
be in a very narrow region as the quasiparticle approx- 
imation fails when r£_F ~ 1. In addition, changing the 
sign of the scatterer changes the sign of the skew effect 
and no scaling is observed. Although this simple model 
seems to capture qualitative aspects of the three regions, 
to make a quantitative link to 3D materials with much 
more complex behavior seems premature at this stage. In 
our calculations, we also identify the hybrid skew scat- 
tering regime of the AHE resulting from the higher order 
scattering processes. Such processes appear to be impor- 
tant for the AMR as the nonvertex and vertex diagram- 



matic parts cancel each other for the Rashba modeli^ 
Our results suggest that the relaxation time approxima- 
tion is not always reliable for the AMR calculations as it 
has been shown recently within the Boltzmann equation 
treatment 1^ 

The paper is organized as follows. In Sec. II, we de- 
velop a general formulation of transport in multiple-band 
non-interacting Fermi systems with further generaliza- 
tions in Appendix A. In Sec. Ill, we calculate the AHE 
in 2DEG ferromagnet with spin-orbit interaction. The 
analytical and numerical results are followed by discus- 
sions and comparison to other works. In Sec. IV, we 
calculate the AMR in 2DEG ferromagnet with spin-orbit 
interaction. Finally in Sec. V, we present our conclu- 
sions. 



II. TRANSPORT IN MULTIPLE-BAND 
SYSTEMS 

The method presented in this section can be applied 
to a multiple-band system described by a Hamiltonian 
Hf) + V{r) that is a matrix in the band (chiral) index. 
In this section, we first derive general non-linear equa- 
tions using nonequilibrium diagrammatic technique, fur- 
ther restricting our consideration to a linear response the- 
ory. 



A. Quantum kinetic equation 

We start by defining the following Green's functions;^ 

Gn = -I (rc*(i+)*t(i+)) = (r*(i+)*t(i;)^ , 
G21 ^ -I (T,vi/(i_)vi/t(i;)^ = _j (vi/(i_)vi/t(i;)^ , 

G12 = -I (T,*(l+)*t(i'_)^ = i (*t(i'_)*(i^)) , 
G22 = -t (Te*(l_)*t(i'_)^ ^ ^T*(l_)*t(i'_)^ ^ 

(1) 

where Tc is the generalized time ordering operator acting 
on the Keldysh contour which can be split in two time 
axis t+ (forward) and t- (backward), ^ is the vector in 
the band (chiral) space corresponding to the Fermi field, 
and 1± = (r, t±) is the variable that describes the spatial 
variable r and the time variable t. The generalized time 
ordering operator performs an ordinary time ordering T 
for the time t+, an anti-time ordering T for the time t_ 
and in the mixed case t_ occurs always after t+ within 
the Keldysh time contour. We can now define the Green's 
function in the Keldysh space: 

G = ( - ) (2) 

\ G21 G22 / 



The scattering potential due to impurities in the Keldysh 
space has the form: 



V(l 1') = ( ^ 



(3) 



where V{y) describes the potential in the band (chiral) 
space formed by many scatterers which for current con- 
sideration can have any general matrix form. The nega- 
tive sign arises here simply because the lower branch in- 
tegration is taken from -l-oo to — oo while in the Keldysh 
loop the time goes from — oo to +00. The Green's func- 
tion in Eq. ^ allows for a perturbation expansion re- 
lying on the Feynman rules. However, the four matrix 
elements of a so defined Green's function are linearly de- 
pendent, i.e. G12 + G21 = Gii+G22- Hence it is advanta- 
geous to perform a linear transformation in the Keldysh 
space to eliminate one matrix element in Eq. ((2]): 



G 



1 

1 -1 



Gil G12 
G21 G22 



1 

-1 1 



G^ G< 
G^ 



which leads to the following scattering potential: 



V 



1 

1 1 



V 



1 

1 -1 



1 
1 



V{v)5{l-l'), 



where & = Gn — G12 is the retarded Green's function, 
G^ = G12 — G22 is the advanced Green's function and 
G^ = G12. There are other choices for the linear trans- 
formation, and our choice is dictated by the fact that the 
Green's function G< can be immediately related to the 
distribution function in the Boltzmann equation4^ 

As of now, it is assumed that ^(r) describes some dis- 
ordered potential and all Green's functions are averaged 
over this disorder. In the transformed Keldysh space, the 
Dyson equation^ becomes: 



Go 







Gq ^ 



G^ G< 
G^ 



= 1, 



(4) 



where R, A, and < respectively stand for the retarded, 
advanced and lesser components of the disorder averaged 
Green's functions and self-energies. The symbol ® de- 
notes a convolution (in position, time and band/spin). 
The diagonal components of Eq. ^ , yield the two equa- 
tions for the retarded and advanced Green's functions: 



(Go ' - S""/^) ® G""/^ = i. 



(5) 



The off-diagonal component of Eq. ^ yields the ki- 
netic equation (sometimes called quantum Boltzmann 
equation) which contains the non-equilibrium informa- 
tion necessary to study transport: 



[G^ 



G" 



■G' 



0. 



(6) 



In order to solve Eq. ([6]), one has to calculate the self 
energy of the particular problem. Here we focus on 



rii rii rii rii 
U ilGU 



+ 



Figure 1: The non-equilibrium self-energy calculated using 
the self-consistent T matrix approximation in Keldysh space. 



scattering by randomly distributed identical impurities 
at zero temperature with 



y(r) -^77C/(r-r,), 



(7) 



where describes the positions of random impurities of 
density rii and fj is some matrix in the band index (e.g. 
in section III, it is a unit matrix corresponding to scalar 
impurities, and in section IV, it is a combination of unit 
and unitary matrices corresponding to charged and mag- 
netic impurities) . A common approximation to this prob- 
lem is the self-consistent T-matrix approximation (TMA) 
which takes into account all the non-crossing scattering 
events from single impurities (see Fig. [T]). We assume 
here that the system is uniform and G depends on the 
difference of spatial variables (r — r^) (however, this re- 
quirement can be lifted for the short-range disorder as it 
is shown in Appendix A). In this case, we can sum up 
the infinite series of diagrams in Fig. [T] arriving at the 
following expression for the self energy in the momentum 
representation (for the sake of compact form we use the 
momentum representation here): 



(k|S|k') ==n,(k|f|k)<5(k-k'), 



(8) 



with the following expression for the T matrix operator 
of impurity placed in the origin: 



f={V + V(g>G(g>V+...), 



(9) 



where V = 



U{r)6{l - 1'). Combining the T- 



fi 
r) 

matrix structure T = V(E)[1 + G(E)T] and solving for the 
off-diagonal component we obtain the equation for the 
lesser component of self energy: 



k|E<|k' 



klT^ 



>G^ 



.f^lk 



S{k~k'). (10) 



The retarded and advanced T-matrices are given by 
the usual form 



= (l + f^(^)<gG«(^))<g-l/. 



(11) 



Equations ^ and l|10p form a general closed set of equa- 
tions for G^. In order to solve these equations, we can 
further simplify them by looking for a solution of the 
form 



G" 



G2< 



G< 



(12) 
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where 



(13) 

and the operator np is the Fermi distribution function. 
In the case of zero temperature, np is the step func- 
tion in the frequency representation np{uj) = 0{—uj) and 
np{t, t') = ij [27r(t — + iO)] in the time representation. 
Equations (fT2|) and (fT3|l will allow us to separate the 
Fermi sea and Fermi surface components of the lesser 
Green's function. By substituting Eq. (fT2)) into Eq. |[6l, 
we obtain the kinetic equation for : 



[G- 



Gf-S<®G^= HQ®np ®G^, (14) 



with 



k|S<|k' 



klT^ 



k') 



Gf ®T'4|k^J(k 

where [■•■'?•■•] stands for a commutator. In order to de- 
rive Eq. lfT4|) . Eqs. ([5]) and (flTI) are used along with the 
fact that f^®G^ = np®f^ -f^® np, and 



5G< 



n;^ (Ki 2." - 2." (K) nj?, (15) 

which is a consequence of Eqs. (|10lll|13p . 

The lesser Green's function contains all the informa- 
tion about the transport properties of our system and 
the charge current density can be calculated as: 

J.(y,z) = |Tr (*t(i')o^(^^^^(l)*(l)^^^^^ +C.C. 



y.iyM^) + yUy,4i ))G<(i,i) 



where ■0(1) 



-ihVi 



eA{l)/c^ /to, 'Ot(l') 



1=1' 
(16) 



^iftV]^' — eA. (1 )/cj /to and A(l) is the generalized 

vector potential matrix in the band index that also de- 
scribes spin-orbit interactions; e — — |e| stands for an 
electron charge. 

B. Linearized Fermi surface contribution 



The kinetic Eq. I|T4|) has not assumed linearity in elec- 
tric field strength nor any particular temporal depen- 
dence. Higher order terms in the impurity density rii 
corresponding to non-crossed diagrams have been taken 
into account as the retarded and advanced Green's func- 
tions in Eq. (fT4|) are calculated self-consistently. In the 
following, we solve the problem for linear response theory 
of a uniform and stationary system in the presence of a 
uniform electric field. 

In the presence of slowly varying perturbations, it is 
useful to perform the Wigner transformation, viz. the 
center-of-mass coordinates {X = (R, T)) and the Fourier 
transform with respect to the relative coordinates (fc = 
{k,uj)). However, the Wigner coordinate k associated 
with the momentum operator — zV is not gauge invari- 
ant and consequently it is not the correct choice for de- 
scribing our system. On the other hand, the kinetic mo- 
mentum k(T) = — iV — eAE{T)/{Hc) is gauge invari- 
ant, and as it will be shown below, for the stationary 



case all time dependence can be conceived in k(T); here 
the vector potential Ae{T) describes the external elec- 
tric field. The time derivative within the canonical co- 
ordinates (marked by wave) becomes a combination of 
time and momentum derivatives within the kinetic coor- 
dinates: df = dr + 9Tk(T)9k, = 9r, = 9k and 
i9ii, = d^. 

In the Wigner representation with the kinetic momen- 
tum, the convolution of two operators is approximated 
as: 



A(g}B 



exp 



AB 



A{X, k)B{X,k) 



xAdkB - dkAdxB 



where we use the four vector notations dxdk = diidk — 
dfd^ and df ^ dr + ^dk- Here, we assume that a 
vector potential Ae{T) = — cET which corresponds to 
a uniform electric field E. The first order gradient ex- 
pansion is sufficient for the Hnear response theory, while 
the second order gradient expansion may be necessary 
for time dependent problems and when the Hamiltonian 
Hq is spatially dependent in order to account for the cor- 
responding Berry curvature effects.— Since we are seek- 
ing homogeneous solutions both in space and time with 
respect to the center-of-mass coordinates, the only sur- 
viving terms in the expansion are 



A(^B» AB- jj-eE [di^Ad^B ~ d^Ad^Bj (17) 

Applying the above Wigner transformation to G^ in 
Eq. (fT3l) . we obtain directly 

G< ^np{G''-G^) + ^d^npeE{di,Gf^ + d^G^A, (18) 



2h 



where G^J"^ are the Green's functions evaluated at equi- 
librium, i.e. E = 0. G2 solves the Kinetic Eq. ([6]) up 
to zeroth order in the electric field E, and therefore the 
expansion in E of Gf and Sf starts from the linear in E 
terms. With this knowledge, we apply the Wigner trans- 
formation to Eq. (fT4| , and find the self-consistent simple 
form of the kinetic equation for Gf : 

G< = G« S<G^, - ieEid^np)G%vG^^ (19) 



/(^7^e^(k,k')Gf(k')r4(k',k) (20) 



where i) = dHo/dKk, and T^q^"^ are self-consistent T- 
matrices evaluated at equilibrium. In the following sec- 
tion, we show how to solve the kinetic Eqs. lfT9|) and l(20|) 
for a simple system described by the Rashba Hamilto- 
nian. Whereas solving Eqs. (fT9| and (|20| require only 
the equilibrium retarded and advance Green's functions 
and T-matrices, note that for G^ we need to solve these 
Green's functions up to linear order in E (see below). 

Prom the equations above, it is natural to decompose 
the contributions to G^ into the Fermi sea and Fermi 
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surface contributions^ such that = Gf + = Gf + 
Gfj where 

= G< + j^^{d^nF)e^{d^Gf^+ d^G%), (21) 

G< = npiG^-G"). (22) 

Next, we Hnearize Eq. lfT6|) in E, carry out the Wigner 
transformation and insert the two components of , ar- 
riving at the two corresponding components of the cur- 
rent density: 

r [ dPli did / \ 

Jxiy.z) = J 2^Tr [Gj j , (23) 

^iiv.) = / (^^Tr (g< ^.(,..)) , (24) 

where the Fermi surface {j^i^y ^j) and Fermi sea 
contributions are identical to ones defined within Kubo- 
Streda formalism 4^ Equations ^ and (HU are the 
main results of this subsection. 



III. AHE IN RASHBA SYSTEMS 

In this section, we apply the above formalism to 2DEG 
with exchange field and spin-orbit interaction. A gen- 
eral numerical procedure is followed by analytical results 
valid in the metallic regime in the limit of small impu- 
rity scattering broadening H/t with respect to the Fermi 
energy ef- We end the section with a discussion of the 
numerical and analytical results comparing them to other 
approaches. For convenience, and in order to keep the ex- 
pressions more concise, we introduce here the dimension- 
less units that can easily be transformed into dimensional 
units by following equations at the beginning of this sec- 
tion. Note that our formalism cannot be used close to 
the energies lu = ±h in Fig. [21 as kpl {I is the mean-free 
path) can become very small and the non-crossing ap- 
proximation in Fig. [T]may fail. Nevertheless, we do not 
expect large corrections to our results around these sin- 
gularities as the non-diagonal conductivity seems not to 
be strongly affected by including the crossed diagrams.— 

A. Calculational procedure 



C. Linearized Fermi sea contribution 

In order to calculate the Fermi sea contribution using 
Eqs. ((22I) and (|24|) . we expand the retarded (advanced) 
Green's function and self-energy up to the first order in 
E following the procedure of Onoda et al.i^ 



qR{A) 



G^(^)+eEG^(^) 



0(i?2), 



(25) 



where Gf = iaEG^b^o, = ^3^^^^=^ and Cf/^ 

{Y.eq^'') are the Green's functions (self-energies) evalu- 
ated at equilibrium, i.e. E = 0. The Fermi sea lesser 
Green's function Gfj calculated up to the first order in 
the electric field E becomes: 



G< - npiGi^ - Gfg) + npenoi - Gf). 



(26) 



We now substitute Eqs. ([Ml into Eqs. ^ and ifTTj) 
only retaining linear terms in E in order to arrive at the 
following self-consistent equations: 



G 



(to) 



G^EeG^^ 



G?„ f ^)+ankSf„ 



d^G^q ( v+dnk 



yR ] fiR 

eq j ^ eq 



(27) 



We restrict ourselves here to 2DEG Rashba Hamilto- 
nian with an exchange field h (breve accent here means 
that h is in dimensional units) in order to obtain simple 
analytical results that connect directly with other micro- 
scopic linear response calculations>i^i^>^ 

Hr = i(nk)V2m + ak-axz-h&^ + iV{r), (29) 

where a is the strength of spin-orbit interaction, & are 
Pauli matrices, ftk = —ifiV — eA/c, A{t) = —cEt de- 
scribes the external electric field and V(r) describes the 
impurities. From symmetry considerations, the most 
general form of the Hamiltonian in Eq. ((29l) should 
treat the coordinate r as an operator r 4- fso(fc) with 
fso(k) = Xa X k originating from the projection proce- 
dure onto the band under consideration.— The spin-orbit 
interaction can also include higher e.g. cubic terms rele- 
vant for the bulk InSb and the HgTe quantum wells with 
an inverted band structure.— 1^ Here, only linear terms 
with Rashba symmetry are considered with fso(k) being 
disregarded as we expect effect of Hso = fso(k)VF(r) 
on the AHE to be small for wide band semiconductors in 
which A is relatively small.— The disorder in the system 
is modeled by impurity delta-scatterers: 



(30) 



t^^^\.)^n^J A^f«(^)(k,k')Gr'(k')T.rnk',k), 

(28) 

where in Eq. ^ we also performed the gradient expan- 
sion. Equations l(26|) . l(27|) and l(28l) are the main results 
of this subsection. 



where f ^ describes the positions of randomly distributed 
impurities of density fii . 

We rewrite the Hamiltonian in dimensionless quanti- 
ties: 



^^iik^ 
Sf 2 



ifk • X z — ha-. 



iVoyjir-r,), (31) 




Figure 2: Electronic band dispersions of the Rashba model; 
throughout the paper, ef is the Fermi energy measured from 
the bottom of the lower band while cj_f is the Fermi energy 
measured from the middle of the gap (region (ii)). 



where e^? is the Fermi energy measured from the mini- 
mum of energy, k = Hq is the dimensionless momentum. 
The dimensionality can be restored by substituting ex- 
pressions for the dimensionless units into the final formu- 
las: 



lo 



I h? 



a 



a 



toVo 



rii = hilff, k ~ klo. 



Also note that whereas ef is measured from the bottom 
of the lower band, in the notation below, we introduce 
ujp which is the Fermi energy measured from the middle 
of the gap (region (ii) in Figl2|. 

In the following, we solve Eqs. (flOl) and l(20|) in order to 
find the non-equilibrium Green's function Gf describing 
processes at the Fermi surface, and Eqs. I|27p and l(28|) 
for the non-equilibrium Green's function - primarily 
Fermi sea contribution. 



We calculate S 



RiA) 
eq 



and the Green's functions G 



R{A) 
eq 



using the self-consistent TMA, i.e. diagonal components 
ofEq. 



fR{A) 
eq 



(32) 



(33) 



G. 



R{A) 



eq 



RiA). 



eq 



akyCTx 



where 7 



R{A) 



.R{A). 
^eqO ) 



[h 



X 2 A) 



^eqZ I 



a 



2V2 



(34) 



(k,c.) 



/ dPk/{2^fGiq 
^R{A)^^^ We calculate self-consistent value of the self en- 
ergy Y,^q^\Lo) for each lo by performing sufficient number 



of iterations in Eq. l(33|) in order to achieve the prescribed 
accuracy (see Appendix B for details). 

With the knowledge of the equihbrium Green's func- 
tion G^(k,a;), we can calculate the local densities of 



R 

E 



+ /' \ + /i\ + /'\\ + 



= f + p + ^ + M + + + 

181 8* ^ «i «i iSi^j*8i b) 

Figure 3: (color online). An infinite set of diagrams represent- 
ing the self-consistent TMA in calculating; a) the retarded 
(advanced) self energy - Sf, and b) the lesser component 
of self energy - in Eqs. l[33|l and ll36ll . respectively. 



states: 



D{uj) 



and the total number of electrons: 



N 



dLoDiuj). 



(35) 



The number of electrons changes as we increase the dis- 
order, and following Eq. I|35p . is always adjusted so 
that the total number of electrons is constant. 



The same TMA is also used to calculate S< and S 
from Eqs. l(20|) and (|28l) . respectively 

(Pk 



RiA) 



T«(c.)G<(k,c.)T4(c.) (36) 



.R{A) 



(Pk 

(2^ 



tr'(-)G^*^^(k,-)f«(^)(..) 



(37) 

The TMA with self-consistent calculation of the equilib- 
rium Green's functions Cfq^^ described in Appendix B 
allows us to take into account higher order non-crossed 
diagrams in the concentration of impurities Ui, with 
weak localization diagrams being disregarded. The pro- 
cedure of calculating the retarded (advanced) and non- 
equilibrium self energies in Eqs. ((33l) and l(36|) is repre- 
sented graphically in Fig. [31 In this graphical representa- 
tion, the bold arrow corresponds to the self-consistently 
calculated retarded (advanced) Green's function. 

For the delta scatterers, T-matrix does not depend on 
momentum k which allows us to perform momentum in- 
tegrations in Eqs. ([SG)) and l(37)) . It is then useful to 
introduce the following 2x2 matrices: 



p{uj) 



.R{A), ^ 



dpk 



Gf (k,w) 



(27r)2^E 



R{A) 



(k,C.) 



(38) 



(39) 



The elements of matrices p and pe satisfy a system of 
linear equations obtained by integrating in momentum 
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space the left and right hand sides of Eqs. I|T9|) and (|27|) . 
respectively: 



(40) 



/ 



<Pk A,R(A)rf,RiA), N--R(A). .rf,R{A), nA-R(A) 



(2^) 



2 J- eg 



2' (27r)2 



«('4)-^H(A)\ 



(41) 

The momentum integrations in the right hand side of 
Eqs. (|40|) and l|4ip are done analytically using the general 

form of the Green's functions Gf'j'^\\i, uj) in Eq. ((34l) . 
Without loss of generality, we take the electric field E 
along the y axis E = {0,Ey) and solve the system of 
linear Eqs. l|40|) and (j4l|) for the elements of matrices p 
and Pe in Appendices C and D, respectively. 

With this, we calculate the current from Eqs. ((23l) and 
(|24)) . respectively, with a use of Eqs. JIS]), {21]), JM]) and 
1271): 



7^ 

■Jxiy) 



-ie f J-L—Ti-id'^f'^df^G'^i 

1 



<:iy) 



(42) 



Jx{y) 



„ (Pk du} 
le J ^9 —ehiUp IT 



(27r)2 27r 



^2 [peq'^^^^eq 



■eq 



eq 



where we use dkGeq^^ 



RiA) ^^R{A} 



G^q'^'vG. 



eq 



+ c.c. 

(43) 

which holds for 
the model of delta impurities. In Eq. (|42l) . we perform 
analytical integrations over momentum k and energy uj 
while in Eq. (|43|) , we only perform analytical integration 
over momentum. The results of these integrations are 
given in Appendices E and F for Eqs. l(42|) and (|43|. 
respectively. 



B. Analytical results in the metallic regime 

In the metallic regime, we are able to obtain analytical 
results as it is sufficient to consider only finite number 
of terms in the expansion with respect to the strength of 
impurity in Fig. [51 For the same reason, we are also able 
to generalize the disorder in Eq. ijSO]) (generalization of 
the theory is given in Appendix A) as follows: 



V{r) 



V^Sir - r,). 



(44) 



where is random, the strength of each impurity has 
the same arbitrary distribution and all strength distribu- 
tions are independent leading to the first four cumulants: 

= 0, n^m'),,,, = V,, n,{{V^f),,^ = Vs and 



dis 



Vi where rii is the concentration of im- 



purities. For the disorder described in Eq. ((30l) . we have 

Vb and for the tele- 
as it is men- 



graph white noise disorder we have Vz 
tioned in Appendix A. 

In this section, we first expand the retarded (advanced) 
self energy in Eq. ijSS]) up to the third order in Vo (or up 
to the terms V3 in Eq. {AH). The lesser component of 
the self-energy in Eq. ((M]) has to be expanded up to the 
fourth order in Vq (or up to the terms in Eq. (|A5P ) 
which corresponds to the four legged diagrams in Fig. 
[Sja) . This ensures that the expansion of the conductivity 
ul.y following from Eqs. I|E3IE5P captures all possible 
terms proportional to 1/Vb and 1. 

The expansion of ol\^ following from Eq. (IF3| IS some- 
what simpler as it only contains the terms proportional 
to 1 and its calculation requires consideration of only 
one bare bubble diagram {e.g. summation of vertices 
leads to higher order corrections). In our discussion, we 
thus concentrate on the diagrams for calculating a\.y and 
also present the result for the bare bubble diagram of 



al,y. Note that in the expansion of cr •', it is important 
to properly consider the branch cut of the "In" function 
taken as (— oo,0]. The diagrams in Fig. [H have direct 
correspondence to the Kubo formahsm diagrams in Fig. 
Hlused in Ref. [13. This allows us to separate the conduc- 
tivity into terms that directly relate to each diagram in 
Fig.H 

We distinguish three regimes for the position of the 
Fermi energy with respect to the gap of the size 2/i; (i) 
u)p > h, (ii) —h < ujp < h and (iii) top < —h (see 
Fig. [2]). To simplify formulas, we introduce the following 
notation: 



A± = ^{ak±y + h'\ Xp = ^2ojpa^ + h^, 



duj{k) 



dk 



ojp > h 



A± ± a2 ' 

— h < LUp < h 



|Ad 



LOp < —h 



where h'± is the density of states at the Fermi level and 
k± are the two Fermi wave numbers for the regimes (i) 
and (iii). In the regime (ii), k+ becomes pure imaginary 
and only k- has the meaning of the Fermi wave number. 
Further, we introduce the following parameter: 
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lo-Q + ^zCTz, with 7 



where 7 = J d'^k/{2TrfG§ 
Y + *7'j 72 = + *7l- Note that the two dimensional 
integral over momentum diverges and 7 is calculated by 
introducing the momentum cutoff, see Appendix B. By 
expanding the result of Appendix B up to the zeroth or- 
der in the strength of impurities, we obtain: 



obtain: 



/(«:) 



2 /ha^iy^ 4hkW Shkia"^ Skta'^X'i 



8h{2h'^ + 2LJFa^ 



kl)Y 



2hkW Vi 



-A 



{kl - 2ujf) In 



fc2 



{kl - 2[j_F)ln 



ftp 



7r(fc2 



In 



4 1 



kl{kl 



7 = ' 



, uop > h 

— —h < ujp < h 



kl + k\- Aujp 



ujp < —h 



2{ki - kl) 

-- — , -h<u:F<h , 
2h 



where k^ is the cutoff in the momentum integration. 

As it follows from the Appendices E and F, the non di- 
agonal conductivities cri^y^^ can be calculated by properly 
choosing the "In" branch that corresponds to the regimes 
(i), (ii) or (iii), respectively. The result of expanding 
Eqs. I|E3|E5P and Eq. ((F3|) for conductivities a^y and 
cr^^, respectively, in the region (i) {up > h) becomes: 



2e^a' 
hi: 



A 



e^ha^ In 


k+{kQ 
kl{kl- 


kl) 
kl) 


hTT'^{kl~kl) 



(45) 



Txy 



0, 



which reproduces result of Ref. |2J in the limit of large 
cutoff fco. In reference to the Kubo formula formalism, 
we can claim the following: the diagrams in Fig. SK) van- 
ish after summation (the intrinsic and side-jump contri- 
butions defined in Ref. \^ cancel each other) ,— the dia- 
grams in Figs. [4j3)-d) are all proportional to t-^ — t — = 

A-i- A— 

and also vanish, and the diagrams in Fig. [4^) lead to the 

result in Eq. ([45]). 0":^^^'^ is zero as the corresponding bare 
bubble contribution in Eq. (|F3|) vanishes. Repeating the 
same procedure for the region (ii) {—h < up < h), we 



ii(ii) 

(7 xy 



-(1 



(46) 

where the diagrams in Fig. [4^) lead to the first three dis- 
order independent terms in Eq. (|46|) (the intrinsic, the 
side-jump and the disorder independent skew scattering 
terms, respectively) pi the skew scattering diagrams in 
Figs. Ilia) and e) lead to the term in Eq. (|46)) propor- 
tional to A, and the diagrams in Figs. S]:) and d) lead to 

the terms in Eq. l(46|) proportional to V^/V2. ^xy is 
calculated from a bare bubble contribution given by Eq. 
(|F3[) and also corresponds to the intrinsic contribution. 
Finally for the region (iii) {cjp < —h), we obtain: 



'xy 



32hujW 



+ 



e 

4^ l^(/i2 + a4)2(V _ ^.2 ) 

MK^ujp + 2a2/i2 - ■iujpa'^) . 

-Iz + 7 



a2(fc2 



{K^ -uopa^){h^ + a^) 



ii{iii) 



32{h^ 
ha^{kl - 



kl)' Vi 



A{h2 + a^)^ 



/i(A_ - A+) 



Ann {a^ - X-){a^ - \+)' 



(47) 

where the diagrams in Fig. [4^) lead to the disorder in- 
dependent term in Eq. (|47l) (it includes the intrinsic, the 
side-jump and the disorder independent skew scattering 
contributions), the skew scattering diagrams in Figs, lib) 
and e) lead to the term in Eq. (|47|) proportional to A, and 
the diagrams in Figs. S]:) and d) lead to the terms in Eq. 

(|47l) proportional to V^/V^- cri^^"^ is again calculated 
from a bare bubble contribution given by Eq. IjFSp . 

The diagonal conductivities can also be calculated by 
expanding Eqs. (|E4|) and fEe]) : 



Up 



h nV2 
h 'nV2vlnl ' 



Up > h 



-h < LUp < h 



{LOp + a^){a^ + Xl) 
h TTV2{a'^ + h^) ' 



LUp < —h 



where we only present the dominant non- vanishing terms 
V^^ as the higher order terms are quite cumbersome. 



C. Numerical results and discussions 

Here, we present results of our numerical calcula- 
tions based on the formalism developed in Section III A. 
Figures I5I6I7I and [8] show the numerical results for 
the anomalous Hall conductivity as a function of the 
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Figure 4: Different diagrammatic contributions to aiy within 
the Kubo formula formalism; a) the ladder diagram (ver- 
tex) contribution ~ 1, b) the skew scattering contribution 
~ l/(n,iVb) ~ V3/V2, c) the double skew scattering contribu- 
tion ~ l/m ~ Vi /Vi, d) the skew scattering contribution in 
which the retarded (advanced) self energy is calculated up to 
the third order ~ l/ui ~ V3 /Vi and e) the fourth order skew 
scattering contribution ~ l/rii ~ V4/V2. 




Figure 5: The anomalous Hall conductivity axy as a function 
of the Fermi energy ujf and the Born scattering amplitude 
7Born. The parameters are chosen as 2c? jh = 35.9, fco = 12 
and Vb = 0.1. The Fermi energy ldf corresponds here to the 
clean system and it is renormalized according to Eq. I|35p in 
the presence of disorder. 



Fermi energy ujp and the first Born scattering amplitude 
iBorn = niVQTTi. The strength of the spin-orbit interac- 
tion is chosen to be the same as in Ref. [U, 2a^/h = 35.9 
(2a^/£'res = 3.59, i^res = Wh; and the strength of im- 
purity is Vb = 0.1, 0.3, —0.1 and —0.3. For the retarded 
(advanced) self-energy, the cutoff in the momentum inte- 
gration is ko = 12 which corresponds to the energy cutofl 
of Ref. m, £c = 3i?res- The Born scattering amplitude is 
varied by changing the impurity concentration rii. 

In the clean limit, when ^Bom 0, we observe skew 
scattering behavior {(Jxy ~ l/riiVb) in which \(Jxy\ rapidly 
increases. For repulsive scatterers (Vb > 0, see Figs. [5] 
and m , the negative conductivity diminishes as we in- 
crease the Fermi energy, until the point lof = —h is 
reached. At this point, the conductivity suddenly in- 
creases without a change of sign, in contrast to Ref. [l| 
where the sign change has been observed but in agree- 
ment with Refs. [13 and [13 (note that Fig. [5] is calcu- 
lated for exactly the same parameters as Fig. 5(c) in 
Ref. [l|). As we increase the Fermi energy further, the 
conductivity increases again around ojf = h acquiring 
a very small negative value. In this regime, both sub- 
bands are partially occupied and only the higher order 
skew scattering22ii21 (hybrid skew scattering) contributes 
to the anomalous Hall effect. Relatively large hybrid 
skew scattering is present in Fig. [6]compared to Fig. [5]as 
the hybrid skew scattering contribution is proportional 
to l/rii ^ Vb^/7Born and should be larger for greater 
impurity strength.— The same is true for the conven- 
tional skew scattering proportional to l/VoUi ^ Vo/jBom 
, which can be immediately seen from Figs. [H[6l[7]and 
m For attractive scatterers (Vb < 0, see Figs. [7| and 
[8]) the sign of the ordinary skew scattering dominating 
in the clean Hmit is opposite to the sign of the ordinary 
skew scattering for the repulsive scatterers. The conduc- 



tivity now increases until we reach the point cop = —h in 
which we observe a sudden drop. One more drop happens 
around the point up = h where the anomalous Hall con- 
ductivity changes sign (see Figs. [7]and[8]). This change of 
sign is consistent with the fact that the higher order (hy- 
brid) skew scattering (prevailing when both subbands are 
partially occupied) does not change its sign as we change 
the sign of disorder.— Comparing Figs. [7]and[8l we again 
see that the hybrid skew scattering is more pronounced 
for larger impurity strength. 

As we increase the disorder by increasing ^Bom, the 
skew scattering becomes less important while the other 
mechanisms, such as intrinsic and side-jump, become 
more important. The intrinsic conductivity only gradu- 
ally decreases with the disorder because the only effect of 
disorder on the intrinsic component comes from broaden- 
ing of Green's functions used in the calculation of the in- 
trinsic component. For repulsive scatterers (Vb > 0), the 
skew scattering has sign opposite to the sign of intrinsic 
and side jump contributions in the region —h < up < h 
(see e.g. Refs. [13 and [13). This explains the sign 
change we observe in Figs. [5l [6] and [TO] in the region 
—h < ujp < /i as we increase jBom (more detailed plots 
are presented in Appendix G). 

The positions of points in which the AHE vanishes can 
be estimated by comparing the Fermi sea intrinsic term 
a^l with the skew scattering term in Eq. l(46|) as those 
two are the major contributions. Physically, the AHE 
vanishes because the intrinsic deflection of electrons be- 
tween the scattering events can be balanced by the skew 
scattering events (in the cross-over region between intrin- 
sic and extrinsic mechanisms). As the former does not 
rely on impurities and the latter does (and changes sign 
with impurities changing sign) , we can have full cancella- 
tion of the two by choosing the proper sign and strength 
of impurities. 
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Figure 6: Identical to Fig. [5] plot but for larger strength of 
impurity Vo = 0.3. 



Figure 8: Identical to Fig. [5] plot but for negative strength of 
impurity Vo = —0.3. 




Figure 7: Identical to Fig. [5] plot but for negative strength of 
impurity Vo = —0.1. 



1. Anticrossings and scalings 

As it can be seen from Figs. [H [6l [7| and [8l the anoma- 
lous Hall effect is resonantly increased around the band 
anticrossing which suggests that for more general band 
structures, the major contribution to the AHE also comes 
from the band anticrossings that happened to be in the 
vicinity of the Fermi level.— This view is well justified for 
the intrinsic AHE in the metallic regime {rep >> 1) as 
it follows from the Thouless-Kohmoto-Nightingale-Nijs 
formulai^ applied to the AHEji This leads to the intrin- 
sic AHE conductivity of the order of e^/ (47rfi) within the 
region (ii) in Fig. [21 The full conductivity that includes 
the intrinsic, side-jump and skew-scattering contribu- 
tions seems to also have the resonant behavior around 
the anticrossing for the Rashba model as it follows from 
our analysis. Whereas our analysis justifies focusing the 
calculations on simplified phenomenological models near 
the anti-crossing locations, we emphasize that it is un- 
likely that these would be characterized universally by 
the Rashba geometry rather than by a combination of 
Rashba and Dresselhauss symmetry. 



In the regime of strong disorder, where this resonant 
behavior is not pronounced, the disorder broadening of 
the Green's functions becomes more dominant and the 
axy has dependence that can no longer be expanded cor- 
rectly in powers of t. This expectation can be easily 
seen from the expressions for axy in our formulation or 
the Kubo formulation, in which 



{a\Vx 

a,l3 



(48) 



where |a) are the exact eigenstates in the presence of dis- 
order and the major contribution for dxy in the dirty limit 
comes from interband matrix elements. When expanding 
things in the momentum basis, the denominator is often 
approximated as {En{k) — En'{k))'^ + {H/tY) while the 
matrix elements are evaluated within the disorder free 
eigenstates. Hence, in the limit of large disorder broad- 
ening, the denominator is simply replaced by (ft/r)^) and 
(Txy ~ (this is different for (Jxx as the contribution from 
interband matrix elements vanishes and Uxx ^ t). This 
of course gives an upper bound for the (Jxy ~ scal- 
ing and in intermediate regimes one would expect 77 to 
be lower than 2. In Figs. [§]and[10l we study the AHE 
calculated in the anticrossing region in order to examine 
in detail the universal anomalous Hall effect regimes that 
could be vaHd for more general band structures. 

We now plot in the logarithmic scale (Jxy as a func- 
tion of (Jxx tuned via rii while all other parameters are 
kept constant. In the clean limit, we recover the skew 
scattering behavior {(Jxy ~ l/riiVo ~ ctxx/Vo) and our 
numerical results (bold line) agree well with the analyt- 
ical results (dashed line) obtained in Sec. IIIB. In the 
moderately dirty limit, we observe the intrinsic-side-jump 
regime (axy = const, this regime is more pronounced for 
smaller Vq) in which the side-jump and intrinsic mecha- 
nisms are dominant. All analytical curves (dashed lines) 
asymptotically reach this regime when axx is very small. 
In the stronger disorder regime, as reported in Ref. [ll, 
the numerical curves have downturn for smaller axx ap- 
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Figure 9: The absolute value of the anomalous Hall conduc- 
tivity \crxy\ versus the conductivity a^x for the spin-orbit in- 
teraction strength 20? jh = 35.9. Dimensionality of quantities 
displayed in this plot is restored. 



Figure 10: Identical to Fig. [9] plot except for the disorder 
which is repulsive here (Vb > 0). Note that the conductivity 
(jjjy changes sign around the cusps. 



preaching the third regime in which a^y ^ with 
</? « 1.6 in Fig. [9l However, a universal scaling can- 
not be claimed since for large and positive strength of 
impurities in Fig. [lO] we only observe the reduction of 
the AHE. 

One should keep in mind that the gradient expansion is 
not fully justified close to the line tef = 1 and our results 
are meaningful only for T£p > 1. Furthermore, since in 
this regime the resonant behavior is strongly diminished, 
in realistic three-dimensional systems, the result could 
be more accurately expressed via the averaged matrix 
elements with some appropriate treatment of the disorder 
broadening. 

Although some experimental works claim to confirm 
the scaling axy ^ 0%^ with around 1.5 . 37, 38, 39, 40, 41 
parison of theory and experiments has to be done with 
care since determining a scaling exponent over a single 
decade is often difficult and has led to many errors in the 
past. For example, in DMS ferromagnets (mentioned in 
Ref. [J to support the scaling hypothesis) the change of 
doping will cause change in the impurity concentration, 
in the magnetization and even in the band structure. 
The theoretical calculations only take into account the 
change in the impurity concentration and further assume 
a Rashba symmetry at the crossing pointsj^i^ Note also 
that within the theoretical treatment, the Hall conduc- 
tivity changes its sign for repulsive impurities (Vb > 0) 
in Fig. [To] which is expected as the skew scattering domi- 
nating in the clean limit has the sign opposite to the sign 
of the intrinsic contribution dominating in the dirty limit 
(see Eq. ijiG])). These types of changes of signs have also 
been observed in experimental systems, e.g. DMS,— and 
of course at that stage scaling is not justified. 



IV. AMR IN RASHBA SYSTEMS 

In the previous section, we showed how the formalism 
developed in Section H can be applied to calculations of 
the anomalous Hall effect in multiple-band systems on 
the example of a Rashba system. In this Section, we per- 
form calculations of the anisotropic magnetoresistance 
(AMR) in 2DEG with the in-plane exchange field, spin- 
orbit interaction and magnetic impurities following the 
same formalism. A general numerical procedure allows 
us to rigorously perform AMR calculations in multiple- 
band anisotropic systems. Within the Boltzmann equa- 
tion approach, such calculations are usually performed 
by using the relaxation time approximation in which the 
transport relaxation time r is calculated from the scat- 
tering amplitudes without fully taking into account the 
asymmetries . This approach was improved in Ref. 
I53I by introducing the perpendicular relaxation time . 
However, in some cases this improvement is yet not suf- 
ficient and Vyborny et al. formulated a procedure for 
finding an exact solution to the Boltzmann equation in 
Ref. l35|. Here we propose an alternative approach for 
AMR calculations in multiple-band anisotropic systems 
to the one proposed in Ref. [H. 

We consider here a 2DEG Rashba Hamiltonian with 
additional in-plane exchange field hx directed along the 
X- axis without any loss of generality : 

Hb, = + ak • o- X z - h^a^ - + V{r), (49) 

where now V^(r) describes the disorder corresponding to 
dilute charged magnetic impurities i^i^'^ 

F(r) = Fo(a<To + 'T,)^<5(r-rO, (50) 
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where describes the positions of random impurities and 
we assume that the magnetic impurities are magnetized 
along the exchange field. The quantity a describes the 
relative strength of the electric part of impurity with re- 
spect to the magnetic part. Note that the AMR is mea- 
sured by changing the direction of electric field E which 
is equivalent to changing the direction of the exchange 
field. 

For the AMR, we only need the diagonal conductivi- 
ties, thus the Fermi sea contribution given by Eq. l(43|) 
vanishes. The AMR can be calculated from Eq. (|42)) and 
we only need to calculate Green's functions at the Fermi 



level. We calculate S 



and Green's functions G 



R{A) 



using the self-consistent TMA: 



fRiA) 
eq 



(51) 



0.2 



01 

< 



0.1 



0.0 




/-"V — Vo=-0.9 
\ -.-Vo=-0.6 




0.0 0.5 1.0 1.5 



a/co 



1/2 



i?(A) 



yR(A) 



(52) 



_ ^R(A).^ 



(c- 



2 



-{h- 



^R(A) 
^eqO 



(h 



egZ - (h 



eqZ 



(c. - f - sfJo^))2 -{h- Sf,(/))2 - a2fc2 + 2h,aky 

(53) 

where 7«(^) = J k / (2^)^ Gfj^\k, lu) = 7^(^)ao -t- 
^x''^^ o'x+^f^'"^'' CTz. We calculate the self-consistent value 
of the self energy S^^^^(w_f) by iterating Eq. l(52|) until 
the prescribed accuracy is reached. 

As soon as we know the T-matrix, we can substitute 
it into Eq. l(40|) and find the matrix p by performing the 
momentum integrations in the r.h.s.. Finally, by substi- 
tuting p into Eq. l(42|) we can calculate the conductiv- 
ity. Note that throughout this section, the angular part 
of the momentum integrations is calculated analytically 
while the radial part is calculated numerically. 

The anisotropic resistance in our system is defined as 
follows: 



AMR 



^yy 



'yy 



and it describes the relative difference in conductivity for 
current fiowing parallel or perpendicular to the magneti- 
zation (represented by the exchange field and/or impurity 
magnetization). 

First, we calculate the anisotropic magnetoresistance 
in Rashba system with in-plane exchange field and non- 
magnetic delta scatterers (see Eq. [30l the magnetic scat- 
terers are absent in this model). Kato et al. found van- 
ishing AMR in the regime (i) (see Fig. [2]) when both 
subbands are partially occupied due to the cancellation of 
the nonvertex and vertex parts in the Kubo formulation. 
In Fig. [m we observe the non-vanishing AMR in the 
regime (i) and this suggests the importance of the higher 



Figure 11: The AMR (100% corresponds to AMR = 1) as a 
function of the dimensionless spin-orbit interaction strength 
a/ y/uJ^ {ujp is counted from the middle of the gap in Fig. 
[2]). The parameters of the model are: Vb = —0.6 and —0.9, 
rii/uiF = 0.01 and hx/u!F = 0.3 . 



order diagrams (such as plotted in Fig. [4^)) not only for 
the AHE but also for the AMR. The AMR effect result- 
ing from the higher order diagrams is more pronounced 
for the larger strength of impurities, similar to the AHE. 
The AMR approaches its maximum around the point at 
which the exchange energy is comparable to the spin orbit 
energy, 2hx ~ a. We note that the non-zero but compar- 
atively weak magnitude of the AMR here in the Rashba 
system is reminiscent of the results in three-dimensional 
DMS ferromagnetsi^ This agrees with physical intuition. 
Under comparison of two mechanisms by which AMR can 
arise - carrier polarization/anisotropy in wavefunctions, 
and impurity polarization / anisotropy in scattering oper- 
ator (see Fig. 1 of Ref. [s^) - the former implies a compe- 
tition between the exchange and spin-orbit terms (in the 
Hamiltonian) resulting in reduced anisotropy strength. 

Consequently, even though observation of the AMR 
effect is deemable in the absence of magnetic scatter- 
ers, we expect much more pronounced effect when the 
magnetic scatterers are present. Our numerical results 
in Fig. [T 2| (p lotted together with the analytical results 
from Ref. 1351 ) confirm this. For the case when the Fermi 
level crosses only one band (region (ii) in Fig. [2]), it was 
found in Ref. \m that AMR = 1/(2 - a^) when \a\ < 1 
and AMR = 1/a^ when \a\ > 1, provided the exchange 
fields are small. For the case when the Fermi level crosses 
two bands (region (i) in Fig. ^ it was found in Ref. Issl 
that AMR = when |a| < 1 and AMR = l/a^ when 
\a\ > 1, in the limit of large Fermi energy (compared 
to the spin-orbit and exchange splitting). We observe a 
perfect agreement between our numerical results and the 
analytical results from Ref. [H. The result in Fig. [T2k ) 
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Figure 12: The AMR as a function of the relative strength a of 
the electric and magnetic parts of impurity potential, by solid 
line we plot analytical results and dots represent numerical 
results; a) Fermi level crosses only one band {lof ~ 0) with 
the following dimensionless parameters Vo — 0.05, a — 1.4, 
m = 0.0015, h:, = 0.0015 and/i = 0.015; b) Fermi level crosses 
both bands Vo = 0.05, a = 0.03, m = 0.002, ha, = 0.002 and 
h = 0.001. 



cannot be reproduced within the common approximate 
approaches22iMiS based on the relaxation time approxi- 
mation as it was pointed out in Ref. [H. The non-physical 
divergence in ayy at the point a = 1 in Fig. [12] is caused 
by the special choice of the scattering potentiali^ As 
soon as the spatial dependences of the electric and mag- 
netic parts cease to be identical (in Eq. l(50l). they corre- 
spond both to delta-scatterers) the divergence of the ayy 
is removed (causing AMR< 1). 



V. CONCLUSIONS 

We have developed a framework for transport cal- 
culations in multiple-band non-interacting Fermi sys- 
tems. By applying this framework to Rashba 2DEG, 
we have resolved some recent discrepancies related to 
the AHE in such systems. The findings of this simple 
2D model have been linked to higher dimensional sys- 
tems arguing that most likely the major contributions to 
the AHE come from the band anti-crossing regions sim- 
ilar to one observed in the Rashba model. Our analyt- 
ical and numerical results reveal the crossover between 
the skew scattering dominated regime in clean systems 
(cxy ~ Vo/jBorn ~ o^xx) and the intrinsic dominated 
regime in moderately dirty systems {(Txy ^ const). In 
dirty systems, we observe the third distinct regime also 
dominated by the intrinsic contribution. In this regime, 
the AHE diminishes in a manner similar to axy ''^ cFxx 
with ip being close to 1.6. This, however, cannot be 



called by scaling as the theory is not meaningful in a 
sufficiently wide range of Gxy and Uxx due to breakdown 
of the quasiparticle approximation when tef ~ 1. For 
the repulsive impurities, we observe that the intrinsic and 
skew anomalous Hall effects have opposite signs. As a re- 
sult, the crossover between those two is also accompanied 
by the change of sign of the AHE. We suggest to engi- 
neer samples with repulsive impurities in order to see this 
change of sign in the AHE. 

We have resolved some discrepancies between the AHE 
results obtained by using the Keldysh, Kubo and Boltz- 
mann approaches by considering the higher order skew 
scattering processes. We have also shown that simi- 
lar higher order processes are also important for the 
AMR when the nonvertex and vertex parts cancel each 
other. We have calculated the AMR in anisotropic sys- 
tems properly taking into account the anisotropy of the 
non-equilibrium distribution function. These calcula- 
tions confirm recent findings on the unreliability of com- 
mon approximate approaches to the Boltzmann equation. 
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Appendix A: GENERALIZATIONS FOR 
SHORT-RANGE DISORDER 

In section IIA, we derive the kinetic equation with the 
self-energy expression that is valid for uniform systems. 
Here, we generaHze this self-energy to non-uniform sys- 
tems in the presence of a short range disorder postulated 
by the following infinite set of correlators: 

{VV)^{V) {V)+V25.,r,, 

{VVV)=Y.{VV) {V)+V^5,,,,,„ 

{WW) = E {{VW) {V) + {VV) {VV)) + V4Sr,r,r,r„ 

(Al) 

where we sum all possible decouplings of the correla- 
tors into a product of two lower order correlators and 

i5rir2r3...rN = Ili=i. .N -1 ^i'^^ ~ ^i+i)- ^ote that usually 
the averaged impurity potential is zero, (V) = Vi = 0. 

After performing the averaging procedure for the 
Green's function, we again arrive at the kinetic Eq. ^ 
with the self-energy given by the following formal expres- 
sion: 
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{Vo + Vo(g>G(g>Vo + ...) k,"^v„, 



(A2) 



where Vq = Vq 



(5(1 — 1') and in the term of n- 



th order proportional to Vq" we replace by Vn which 
ensures that the correlators in Eq. IjAip are properly 
considered. It is convenient to introduce the notation: 



E = {Vo + Vo(E>G(g>Vo + ...) , 



(A3) 



which, in analogy with the self-energy, also has retarded 
(advanced) and lesser components. Eq. (jAsp 

can be rewritten in the form of T-matrix equation, E = 
V (E) [1 + G ® E], which leads to the expressions for the 
self-energies: 



{E^'®G<®E^) |y.^y„ 



(A4) 



where the notation |yn_v„ is formal and it means that 
}t;R{A) j^g^g expanded with respect to Vq and 

then the substitution has to be appHed. Eqs. lfT4|) and 
(fT5|) can now be rederived for non-uniform systems with 
the disorder given byEq. (fATI) . 

Nevertheless, for the purposes of this paper, it is suf- 
ficient to consider the uniform and stationary case. This 
leads to substantial simplifications outHned in sections 
IIB and IIC. Results of sections IIB and IIC also hold 
for the disorder given by Eq. (|Aip with the exception 
of Eqs. ((20|) and (|28l) that should be replaced by the 
following equations: 



rG<(k') < 



E, 



R{A) 
eq 



(2^)2 



R[A) 



(k') E 



rR{A) 

-'eq 



(A5) 

Finally, we would like to present several examples in 
which the disorder given by Eq. (jAip is reaHzed. The 
simplest example is given by Eq. (|30)) and in this case 
Vn = Vq . For the disorder given by Eq. (|44| . we have 
Vn = '^j ((K))")tiis- ^^''^ the Gaussian white-noise dis- 
order, only V2 is non-zero and Vn\n^2 = 0. For the 
telegraph white-noise disorder all odd correlators vanish, 

V2n+1 = 0. 



Appendix B: CALCULATION OF 
SELF-CONSISTENT SELF ENERGY t 



The following relations can be calculated by a direct 
analytical integration of Eq. l(34)) : 



{K+ - 2W) [ln(A-o -K+)- H-K+)] 

211 {K^ ~ K+) 
(A'_ - 2W) [InjKo -K_)- \n{-K_)] 

211 {K^ - K+) 



ln(/^o - K+) - H-K+) - MKq - K_) + \n{-K_ 
7r(A'_ - K+)/H 



7 = (7 ) , Ixiv^z) = hx{y,^)) ^v(^) = 



where W ^ uj -Ysfq^. H = h - Ef^^ , K± = 2{W + T 

^JH"^ + 2Wa'^ + a^) and Kq ~ describes the cut-off 
ko in momentum integration. 

For each energy, Sf^g(a;) and S^^l"^) are calculated by 
performing a number of iterations with the consequent 
iteration according to 



n,Vo{i-Voj''^^T'^o 



^R{A) 
-'eqZ 



-Tr 



i,l^o(i-K,7«^^))-i<7, 



The iterations are performed until the prescribed accu- 
racy is reached. 



Appendix C: CALCULATION OF THE MATRIX 

p{uj) 



For the electric field E along the y axis E = (0, Sy), 
we solve here the Hnear Eq. l(40|) for the elements of the 
matrix piuj) by performing analytically the momentum 

integrations of the Green's functions G^^"^'' (k, u) (given 
by Eq. ((34|) ) in the right hand side. For each energy uj, 
we obtain the following expressions that also depend on 
the self-consistent values of E^g(w) and Y^q^{uj): 



p—{uj) = p++{uj) = 0, 
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{uj) ~ id^npeEya 



Kl+A{H + W){H*-W*) ln(-/f_) K*J + A{H + W){H* - W*) \n{-K*_) 



(A'_ - K*_){K^ - K+){K^ - KI)t: (A'_ - K*_){K*_ - K+){K*_ 

4(i7 + W){H* - W*)] \n{-K+) [Kf + A{H + W){H* - W*)] ln(-A:;) ' 



-1 

{in - 



{K^ - K+){K*, - K+){K+ - K^TT (A-_ - A- )(A- - Kl)iK+ - K^-k 

[2H - A'_ + 2W){2H* +K^~ 2W*) lii(-A'_) {2H - A'* + 2W){2H* + K*_ - 2W*) \n{-K*_) 



(A:_ - K*_){K^ - K+){K^ - A:*)7r 
■ 2W){2H* + K+- 2W*) \n{-K+) {2H - K'_ 



(A:_ - K*_)(K*__ - K+){K*_ - a:* )7r 
2W){2H* +KI- 2W*)ln(~Kl 



(X_ - K+)iK*_ - K+){K+ - Kl)n 



[K. - KX){K*_ - KX){K+ - K^TT 



n,T++Tl 



Kl +A{H -W)iH* + W*) ln(-A'_) Kl^ + A{H - W){H* + W*) ln(-A:*) 

'9^^^^Ey^[ \k^-K*_){K^-K^){K. 

{Kl + A{H - W){H* + W*)) ln(-A + ) 
^ (A_ - A+)(A'* - A:+)(A:+ - Kl)-K (A'_ - Kl){K*_ - Kl){K+ - KD-k 

{2H + A"_ - 2W){2H* - A:_ + 2W*) ln(-A'_) {2H + AT* - 2W)i2H* - AT* + 2W*) ln(-A:* ) 



- KI)tt (A:_ - K*_){K*_ - K+){K_ 

{Kf + 4(i7 - W){H* + VK*)) ln(-A:. 




(2i7 + 



(A'_ - K*_){K_ - K+){K_ - A"^ 
2VK)(2i7* -K+ + 2W*) ln(-A'+) 



)7r (AT. - K*_){K*_ - K+){K*_ - Kl)n 

{2H + a:; - 2W){2H* -KI + 2W*)Iti{-KI' ' 



- K+)iKl - K+){K+ - Kl)n 



(A:_ - KX){K*_ - KX){K^ 
I 



Appendix D: CALCULATION OF THE MATRIX 

For the electric field E along the y axis E = (0, Ey), we 
solve here the linear Eq. l(4T|) for the elements of the ma- 
trix Pe^"^' (uj) by performing analytically the momentum 
integrations of the Green's functions Gq ^^\]s., u) (given 



J 



by Eq. ((34|) ) in the right hand side. For each energy uj, 
we obtain the following expressions that also depend on 
the self-consistent values of T,qq{uj) and T,q^{uj): 



0, 



p1,{uj) = Ata H{-1 + a.S^o) - Wd^^g 



-ln(-A'_) + ln(-A"_ 



Kl - K'l + 2K^K 

(A'_ - K+) KlK+Ti + A'3 A'_7r + K'iK+n{r^^T++ - AK+n,T^^T++ {H^ - W^) 
KlK^n,T^^T++ - AK^n,T^^T++ {H^ - W^) - 2nKlKl - 8K+K^n^T^^T++W 
+2K_K+niT__T++ \aH^ + (A'_ - 2VF)(-A:+ + 2W) \ [ln(-A:_) - ln(-A:+) 



Appendix E: CALCULATION OF THE FERMI 
SURFACE CONDUCTIVITY 



For the electric field E along the y axis E = (0, Ay), 
we perform momentum k and frequency w integrations in 
Eq. l(42|) . It is convenient to divide the resultant conduc- 
tivity into two parts; the bare bubble part cr^^^yy) that 
corresponds to calculating only the second line in Eq. 
(|42|) effectively assuming that p{luf) = 0, and self consis- 



tent part (jly'(^yy^ that corresponds to calculating the first 
line in Eq. (|42l) that takes into account correction due to 
self-consistent calculation of p{u>f)- 



'xy 



Isc 
xy ' 



(El) 



(E2) 



yy yy ' yy ^ 
We arrive at analytical expressions for the bare bub- 
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ble contributions to the conductivities a^y and cr^^ that 'Sq-^(ujf) at the Fermi surface: 
depend on the self-consistent values of T,qq{lup) and 

I 



ail _ 2*a^ [-H*{K_ + 2W) + H{K_ + 2W*)] \n{-K_) ^ 2ia^ [H*{K*_ + 2W) - H{K*_ + 2W*)] ln(-A'* ) 
~ (K^ - K*_){K^ - K+){K_ - KX)-K^ ^ 

2ia^ [-H*{K+ + 2W) + H{K+ + 2W*)] M-K+) 
{K_ - K+){K*_ - K+){K+ - KDn^ 



+ 



{K^ - K*_)(K*_ - K+){K*_ - K^TT^ 
2ia^ [H*{Kl + 2W) - H{Kl + 2W*)] ln(-i^;) 



{K_ - Kl){K*_ - Kl){K+ - A'* )7r2 



(E3) 



all _ Kl{K^ -■iK+)Kl\TL{-K^) Kl{K^{K^+K*_)-{K_+iK*_)K+ + 2Kl)\n{-K^) 
^V?^ " ^ 4(A:_ - K+)3{K_ - X* )7r2 4(A:_ - K*_){K_ - K+Y{K_ - if* )7r2 

K_K*J{K*_ - 3KI) \n{-K*_) K*J {K*_ {K*_ + K+) - {K*_ + 3A'+)A'; + 2Kf) \ii{-K*_) 
^ 4{K_ - ATI) (AT* - Kl)H^ ^ 4(A:_ - AT* )(A:* - K+){K*_ - KD^n^ 

AT^ (SA-. - K+)Kl ln(-A'+) Kl {2Kl + K+{K*_ + K+) - A'_(3A"1 + K+)) \n{-K+) 

~ 4(a:_ - a:+)3(a:+ - a:* )7r2 4(is:i - k+){k_ - K+f{K+ - ki)tt^ 

K^Kf{-iK*_ - Kl) ln(-A';) Kf {2K*^ + Kl{K+ + K*^) - Kl {iK+ + KD) ln(-A';) 

^ 4(a:_ - a:;) (a:* - A:;)37r2 + 4(is:_ - ki){k*_ - Kif{K+ - ki)tt^ 

a:_a:+ (a:*2 + a:;2) _ a:2 (a:i2 - atia:; + Kf) - Kl {ki^ - kik^ + Kf) 

2{K- - K+Y{K*_ - KlfTT-^ 
2 {mW*a^ + 2HH* {K^ - 2a^) - Kl{W + W* + a^) + 2K^ {WW* + (W + W*)a^)) lii{-K^) 

2{K^ - K*_){K^ - K+){K^ - A:;)7r2 
2 {AWW*a^ + 2HH* (AT* - 2a2) - X*2 (^w + W* + a"^) + 2K*_ [WW* + [W + W*)a^)) \n{-K*_) 
^ 2{K^ - K*_){K*_ - K+){K*_ - A:;)7r2 

2 {mW*a^ + 2HH* (A+ - 2^2) -~ Kl{W + W* + a^) + 2A% {WW* + [W + W*)a'^)) \n{-K+) 

2(A:_ - K+){K1 - K+){K+ - KI)tt^ 
2 {AWW*a^ + 2HH* (A'; - 2a2) _ Kf {W + W* + a^) + 2A'* [wW* + [W + W*)a'^)) ln(-A:;) 
^ 2(K_ - KX){-K1 + Kl){-K+ + if;)7r2 

2 (-2i72 + [K^ +K+~ 2W)W) {{K^ + K+) (4ij2 + A:_A'+) - m^K+W - A{K^ + K+)W'^) o? 

2(A:_ - A:+)27r2 2A_(A:_ -K+YK+tt^ 

-[ln(-if_)-ln(-/f+)] 

2 (VF(-2A'_A+ + (A_ + K+)W) + {-K_K+ + 2(A"_ + A"+)M^ + AW^) + {K_ + K+ - Aa^)) 

^ 2{K- - K+fir'^ 

2 (-2iJ*2 + (AT* + a:; - 2W*)W*) {{Kl + Kl) {4H*^ + KlKf) - SAT* AT^VK* - 4(A:1 + A:;)W^*2) 

2(A:* - A:;)27r2 2K1{K1 - KD^Kln^ 

-[ln(-A-*)-ln(-if;)] 

2 (W^*(-2A:* AT* + (AT* + Kl)W*) + (-A'* AT* + 2(A'* + AT* + 4M^*2^ ^2 _^ + AT* - 4^2)) 

^ 2(a:* - A:;)37r2 

, m 



where in this Appendix all parameter are taken at the 
Fermi surface: W ^ uip — Soo('^-f)) H = h — Y.q-^{llif) , 
K± = 2{W + T VH^ + 2Wa^ + a'^). 



The analytical expressions for the self-consistent con- 
tributions to the conductivities ai^!^ and alf'^ become: 
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ait; _ an, {p+^T++T*__ {KI + A{H + W){H* - W*)) - p-+T__Tl+ {K^_ + A{H ~ W){H* + W*))) ln(-if_) 



aui ( 


[P+- 








2(/f_ - K*_){K^ - K+){K^ - KDn^eEy 


- W){H* H 


hM^*)))ln(-A'*) 




mil 
H 


KP+- 


-T+- 






2{K_ - K*_){K*_ - K+){K*_ - Kl)TT^eEy 
+ 4(/f + - - p-+r TJ^ (A^ + 4(i/ - 


- W){H* + 


W*)))\tl{-K+) 




aiii ( 


[P+- 


-T+, 


_T*_ 


[Kf 


2(A'_ - A'+)(A-* - K+){K+ - Kl)n^eEy 
+ A{H + W){H* - W*)) - p-+T__Tl^ {Kf + 4(ff - 


- W){H* H 


VW*)))\n{-Kl) 














2(A'_ - A';)(-A-* + Kl){-K+ + Kl)TT^eEy 






(E5) 


ayy" iarii 


{P+- 


-T++ 


Tl_{ 


Kl + A{H + - W*)) + p-+T__Tl^ {Kl + A{H - W){H* + M^*))) ln(-A'_) 


e?lh 

iarii 




-T+- 






2(A:_ - K*_){K^ - K+){K_ - K^TT^eEy 
+ 4(iJ + W){H* - W*)) + p-+T^-Tl^ (A*2 + A{H 


- W){H* 


+ W^*)))ln(-A'*; 


1 


ian. 




-T+ 


_+r* 


-{Kl 


2(A:_ - K*_){K*_ - K+){K*_ - A"* )7r2e£;j, 
+ 4(iJ + - M^*)) + p^+T._TIj^ {Kl + A{H - 


- W){H* - 


^W*)))H-K+) 




+ 

iarii 




-T++T1_ 




2(A:_ - K+){K*_ - K+){K+ - K^n^eEy 
' + 4(iJ + - W*)) + p-+T—TX^ {Kf + 4{H 


- W){H* 


+ W*)))\xi{-Kl] 


1 












2(A:_ - A';)(-A:* + Kl){-K+ + KDn'^eEy 






(E6) 



where again all parameters are calculated at the Fermi the following expressions for conductivities and ayyi 
surface. 



Appendix F: CALCULATION OF THE FERMI 
SEA CONDUCTIVITY 



'^ii - 0, (F2) 



For the electric field E along the y axis E=(0,Aj,),we "yy 
perform momentum integrations in Eq. I|43|) arriving at 



ail"- f Uta^H{l-d^i:§o)(K-+K+)+d^J:g{K+W + K^{K+ + W))] 

J I K.K+{K.-K+Y^^ 

2^2 [4(1 - )i7 + d^J:§,{K- +K+ + AW)] (ln(-if_) - ln(-if+)) 



+ c.c. 



(F3) 



^IIsc 



dun, 



p| )T__T++ (4if2 + K^K+ ~ AW) {Kl - Kl + 2K^K+{- \n(-K^) + H-K+))) 



2K^K+{K. 
+C.C. = 



K+)H^ 



r 



The fact that cr^^"'^ = follows from the identity AH^ 



K^K+ - AW"^ 



0. 



As one can see, cr^^ and cr^y*'^ contributions to the 
Fermi sea Hall conductivity vanish and the non- vanishing 



contribution a^^y depends on the self-consistent values of 
T,qq{uj) and Tjq^[uj) and its calculation from Eq. I|F3|) re- 
quires numerical integration over uj. 



Appendix G: DETAILED RESULTS FOR THE 
HALL CONDUCTIVITY 



components of the AHE conductivity, particularly the 
Fermi sea contribution a^y, the bare bubble contribution 



In order to gain more insight into the behavior of the 
anomalous Hall effect, in Figs. [13] and [14] we plot different 
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Figure 13: The anomalous Hall conductivity crj°* = a^y + aiy'^ + a^y and its components {a^y, (t^J'^, cr^y) versus the averaged 
relaxation rate 1/r — 2ImEoo (defined in Appendix B). The spin-orbit interaction strength is 2ma^/i?res = 3.59 (-Eres ~ Wh); 
the strength of impurities: Vb = 0.01, 0.1, 0.2, 0.3; the Fermi energy ep/E-res = 0.9 for luf = 0, ep/Eres = 0.5 for uip = — 4/i 
and £f/ Exes ~ 1-5 for lof = Qh. Dimensionality of quantities displayed in this plot is restored. 
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Figure 14: Identical to Fig. [13] plot with attractive disorder (Vo = 



-0.01, -0.1, -0.2, -0.3). 
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fj^y (this corresponds to f^Jy''^^ in Ref. [l|) and the self 

consistent contribution cr^^"^ (this corresponds to cr^y ^ in 
Ref. HI). In Fig. [131 we take the same parameters as in 
Figs. 7 and 8 of Ref. HI and we find disagreement with 
Ref. [H in the results for the contribution f^^^ i'^iy'^)- 
The contributions (jlt and a^i perfectly agree with Ref. 

H. 

In the clean limit r — > oo, we see that cr^^^ and thus 
the total Hall conductivity (tJ°* diverge. This diver- 
gence (o'i'y'^ ~ i/niVo in the regions (ii) and (iii) and 
a^y"^ ^ l/rii in the region (i), see Fig. [2| is due to 
the skew scattering. The conductivity cr^y^ also con- 
tains the side-jump contribution which can be best seen 
in Fig. [I3k ) in the sharp peak in the conductivity for 
small l/r. The skew scattering contribution decays much 
faster compared to the side-jump and intrinsic mecha- 
nisms as we go to larger 1/t. As a result, we can expect 



a cross-over between the region dominated by the skew 
scattering and the region dominated by the side-jump- 
intrinsic mechanisms. When both subbands are partially 
occupied (see Figs. [T3f ) and[T4]F)), the higher order skew 
scattering is still present. However, we do not expect 
a well pronounced cross-over as the intrinsic contribu- 
tion cancels the side-jump contribution in the metaUic 
regime (see Eq. (|45|) ). By comparing Figs. [T3f) and 
fUT). one can see that the higher order skew scattering 
(hybrid skew scattering)^ does not change sign when we 
change the sign of impurities. 

When the side-jump-intrinsic and the skew scattering 
components have opposite signs, as in Fig. (fT3| . we ob- 
serve the AHE sign change instead of the cross-over. In 
Figs. fT3k)-d). the skew scattering is negative in the clean 
limit while the side-jump-intrinsic part is positive. This 
inevitably leads to the sign change of the conductivity 
axv as we increase the disorder. 
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